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Abstract 

We investigate the dynamics of solitons in generalized Klein-Gordon equations in the 
presence of nonlinear damping and spatiotemporal perturbations. We will present 
different mechanisms for soliton explosions. We show (both analytically and numer- 
ically) that some space-dependent perturbations or nonlinear damping can make 
the soliton internal mode unstable leading to soliton explosion. We will show that, 
in some cases, while some conditions are satisfied, the soliton explodes becoming a 
permanent, extremely complex, spatiotemporal dynamics. We believe these mecha- 
nisms can explain some of the phenomena that recently have been reported to occur 
in excitable media. We present a method for controlling soliton explosions. 
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Solitons are used in many important technological applications. Among these 
applications we can mention long distance communication systems and soliton 
oscillators in superconducting devices [1]. 

However, under certain conditions, solitons can become unstable [2,3,4]. Such 
instabilities have been called both soliton breakup and soliton explosions. For 
instance, Milchev and coworkers have studied the Frenkel-Kontorova model 
with anharmonic interatomic interactions [2]. They have found that, in the 
Frenkel-Kontorova model with nonconvex interactions between closest neigh- 
bors, a breakup of the kink takes place when the effective amplitude of the 
sinusoidal substrate potential exceeds a certain critical value. 



Corresponding author. Fax: +58-212-5041148; e-mail: jorge@pion.ivic.ve 



Preprint submitted to Elsevier Science 



8 February 2008 



An extensive discussion of soliton dynamics in the framework of the Frenkel- 
Kontorova model can be found in the recent book [5]. 

On the other hand, in Refs. [3] it was predicted that the soliton internal mode 
can become unstable leading to soliton explosions. 

Such instabilities as soliton explosions can affect all the mentioned applica- 
tions. So it is very important to understand all the possible mechanisms of 
soliton explosions in order to avoid them. 

In the present letter we investigate generalized Klein-Gordon equations as the 
following 

<Pu + R (&) - (j) xx -G ((f)) = F(x, t), (1) 

where G ((f)) = —dU (<p) /d(f>, U ((f)) is a potential function with at least two 
minima <f>i,(f>3 and a maximum 2 , such that U (0i) = U (<j)z) = , R(<f>t) 
are dissipative terms, and F(x,t) represents external perturbations. We are 
interested in kinks, that is, topological solitons between the points (f>i and (f>z . 
The famous sine-Gordon and 4 -systems are particular cases of Eq. (1). 

The topological solitons studied in the present letter possess important appli- 
cations in condensed matter physics - they describe domain walls in ferromag- 
nets and ferroelectric materials, dislocations in crystals, charge- density waves, 
interphase boundaries in metal alloys, fluxons in long Josephson junctions and 
Josephson transmission lines, etc. [6]. 

We will present different mechanisms for soliton explosions. We will show that 
in some cases, while some conditions hold, the soliton explosion is permanent. 

A soliton destruction is observed when inhomogeneous space-dependent per- 
turbations are present: 

(f)u + l<f>t-<f> xx -G((f)) = F(x). (2) 

We should say that the zeroes of F (x) are candidates for equilibrium positions 
for the soliton [3]. If F (x) possesses only one zero Xq (F (xq) = 0), then it is a 
stable position for the soliton if (dF/dx) x * > 0. Otherwise, the position is an 
unstable equilibrium. The opposite is true for the antisoliton. The center of 
mass of a soliton can make oscillations around a stable zero of F (x), and it can 
move away from an unstable one. However, when F (x) has an unstable zero at 
X — X Q and additionally, the following conditions holds lim^-too F (x) = =F-Foo, 
Fqo > and Foo is larger than some critical value, then the soliton can be 
destroyed and an antisoliton is "created" in the equilibrium position 
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When the soliton is in an unstable equilibrium position, it is "stretched" by 
the pair of forces that are acting on its body in opposite directions. And there 
is a limit for the magnitude of the pair of forces that the soliton can resist. 
Nevertheless, there is a more subtle mechanism for soliton destruction. When 
a soliton is close to an unstable equilibrium position, many internal shape 
modes of the soliton can be excited [7,3]. If (d 2 F/dx 2 ) x * is larger than some 
critical value, then the first internal shape mode can become unstable and 
this instability can lead to the soliton destruction. In this phenomenon the 
soliton decays into an antisoliton and two solitons. What is interesting in this 
situation is that F (x) can be a localized perturbation. 

Suppose we are interested in the stability of a soliton situated in equilibrium 
positions created by the inhomogeneous force F (x) . Using an inverse-problem 
method [3,8] we construct an exact solution <p k (x) with the topological prop- 
erties of a kink-soliton. Then we investigate the stability of the solution solving 
the spectral problem 

Lf(x)=Tf(x), (3) 

where L = -d xx - [dG (0) /dfl^, T = - (A 2 + 7 A). 
Let us discuss some examples in detail. 
The force 

F (x) — a tanh ((3x) \S + e sech 2 {(3x)\ , (4) 



can sustain a kink-soliton equilibrated at point x — 0. 

When we solve the stability problem (3) we obtain the following eigenval- 
ues of the discrete spectrum (for simplicity we assume 5 = (a 2 — 1) /2, e = 
(A(3 2 - a 2 ) /2 and G (0) = (0 - 3 ) /2: Y n = (3 2 (A + 2An - n 2 ) - \, where 
A (A + 1) = 3a 2 /2f3 2 , n < A. 

The translational mode of the soliton is stable if 2/3 2 A > 1. If this condition is 
not satisfied, this just means that the soliton center of mass will move away 
from the unstable equilibrium position x — 0. This does not necessarily lead 
to the soliton destruction. 

However, if the following condition holds: 

2(3 2 (3A- 1) < 1, (5) 



the soliton first shape mode is unstable. In this case, the soliton can be de- 
stroyed! 
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In the very special (but also very illustrative) case e = 0, we have that for 
4[3 2 > 1, the translational mode is stable. This means the equilibrium position 
x = is stable for the soliton. The soliton center of mass can perform oscilla- 
tions around point x — 0. If Af3 2 < 1 (but 10/3 2 > 1), the translational mode 
is unstable. In this case, the soliton can move away from point x — 0, but 
it conserves its very characteristic shape because the internal modes are still 
stable. However, if 10/3 2 < 1, then the first internal (shape) mode becomes 
unstable. In this situation the soliton can explode. If we apply a spatiotem- 
poral perturbation that periodically (in time) creates this instability in the 
place where the soliton is situated at that instant, then the result will be a 
highly nonstationary spatiotemporal state where the soliton is not allowed to 
recover its original shape. The soliton is in a permanent explosion. 

Another very important example of equation (1) is the sine-Gordon equation 
(i.e. G(4>) = -sin0). Suppose F(x) = 2 {f3 2 - 1) sinh (fix) cosh" 2 (fix). This 
perturbation creates an equilibrium position for the sine-Gordon soliton at 
point x = 0. When we solve the eigenvalue problem (3) for the sine-Gordon 
soliton in the presence of this external force we get the following discrete 
spectrum: T n = [3 2 (A + 2An - n 2 ) - 1, where A (A + 1) = 2/f3 2 . The integer 
part of A yields the number of eigenvalues in the discrete spectrum. For [3 2 > 1, 
the translational mode is stable and there are no internal modes. If (1/3) < 
f3 2 < 1, then the translational mode is unstable (but still there are no internal 
modes). When (1/6) < f3 2 < (1/3) there appears an internal shape mode, 
which is stable. For [3 2 < 2/ [A* (A* + 1)], where A* = (5 + ^/Tfj /2, the first 
internal shape mode becomes unstable. This perturbation can destroy the 
sine-Gordon soliton. 

We should say that a soliton, moving in a medium that is homogeneous ev- 
erywhere except for a zone where the conditions for the instabilities hold, can 
undergo dramatic transient changes. But when the soliton leaves the men- 
tioned zone, it will return to its original steady state shape. 

How can we produce a permanent soliton explosion? We can use time-dependent 
perturbations: 

<h + 7& 'Pxx G (0) = F 1 (x) + F 2 (x, t) , (6) 

where Fi (x) is a perturbation that creates a potential well for the soliton (i.e. 
F (x) possesses a zero x$ such that (OF/ dx) x , > 0) and F 2 (x, t) is a space-time 
perturbation that periodically generates the instabilities conditions. Figure 1 
shows an example of those highly complex spatiotemporal behaviors. In all 
the figures G (<p) = (0 — 3 ) /2. However, similar results are obtained with the 
sine-Gordon and other Generalized Klein-Gordon equations. 

Can we produce permanent soliton explosions without time-dependent exter- 
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Fig. 1. Initial steps of a permanent soliton explosion sustained by a spatiotem- 
poral perturbation as in Eq. (6). Here F± (x) = (1/2) A (A 2 — l) tanh (Bx), 
F 2 (x,t) = (1/2) f A(AB 2 - A 2 ) cos (ujt)smh(Bx) cosh" 3 (Bx), A = 1.5, B = 0.2, 
7 = 0.1, u = 1, / = 3. 

nal perturbations? 

Here we would like to remark that solitons can move with constant velocity 
(without attenuation) in active and excitable media, and in systems with 
nonlinear damping even without explicit external forces. 

Let us discuss here briefly the importance of nonlinear damping. Linear dissi- 
pative systems like the damped harmonic oscillator tt + 70 t + u; 2 = cannot 
sustain oscillations. However, the nonlinear oscillator 4> a — 60 t + a0 3 + a;o0 = 
supports a stable limit cycle [9]. The transition from a stable focus to an un- 
stable focus and a stable limit cycle is the result of a Hopf bifurcation. This 
system is very easy to realize in practice using negative-resistance electronic 
elements [10]. 

Soliton systems as the following: 



where dR(<p t ) /d(p t is negative for small values of \(f> t \ and positive elsewhere, 
can support solitons moving with a constant velocity. An example of this kind 
of systems can be realized in practice using a Josephson junction transmission 
line where the resistor is a negative-resistance twin-tunnel-diode circuit or a 
twin-transistor system [10]. In this case, R (<f) t ) = —b(f> t + a0 3 is a good model. 

We have investigated the shape mode stability of the soliton in the presence 
of nonlinear damping as we did before using the spectral problem (3). 

Suppose R ((j) t ) possesses two local extrema: a maximum and a minimum such 
that the value of \R((j) t )\ at these extrema is R m . If this value is comparable 



<t> tt + R (4>t) ~ 0: 



XX 



G (ct>) = 0, 



(7) 
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Fig. 2. (a) Soliton explosion due to nonlinear damping. Here in Eq. (7): 
R(4>t) = —b(pt + a = 1, b = 0.7. (b) Limit cycle produced with the dy- 
namics of the soliton center of mass in Eq. (8). Here T (x) = 1 — 1/ cosh(Dx), 
F (x) = ^tanh (Bx), A = 0.45, B = 0.65, I = 2, D = 0.65. 

with the absolute value of the extrema of G (</>) (let us call it G m ), then the 
soliton can be destroyed. In fact, if R m > G m the internal shape mode of the 
soliton can be unstable and the soliton becomes a highly nonstationary state. 

When F (x) in Eq. (2) has a stable zero, say x = Xq, the center of mass of a 
soliton can perform damped oscillations around x* . If we wish to sustain these 
oscillations without explicit time-periodic external forces, then we should re- 
sort again to negative damping. Another way to experiment negative damping 
is when the damping coefficient in Eq. (2) is a function of x: 

<P tt + Y{x)<p t - ( t )xx -G{ ( j ) ) = F{x). (8) 



Here V (x) is negative in a neighborhood of Xq and positive elsewhere. This can 
be done in a chain of nonlinear oscillators using negative-resistance circuits 
[10] only in some small interval of the chain. An example of V (x) with the 
required features is T (x) = 7 1 — L/ cosh 2 (Dx) , where (1 — L) < 0. 

Figure 2 (b) shows a limit cycle which is the result of the dynamics of the 
soliton center of mass in Eq. (8). However if we are not careful, the soliton 
can explode also in this system. We have solved the soliton stability problem 
for this equation. The most important result is that the first internal shape 
mode is unstable for L > 5/2. This behavior can be observed in Fig. 3 (a). 

Can we control all these types of explosive dynamics? Well, if we can change 
the parameters of the system, then we can use parameter values that do not 
lead to soliton explosions. However, sometimes we can not change the param- 
eters. We just are allowed to apply some external perturbation. 
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Fig. 3. (a) Highly nonstationary spatiotemporal dynamics produced by Eq. (8). Here 
the simulated equation is the same as in Fig. 2 (b) with / = 6. (b) The soliton dynam- 
ics can be controlled in order to avoid the soliton explosion. Here in Eq. (9) 7 = 0.1, 
F p (x,t) = -0.385 tanh (Bz) cos (wi), F c (x,t) = -0.75 tanh (Bx) cosh~ 2 (Bx), 
B = l,u = 0.2. 

Let us pose the following problem: 

<t>tt + l<t>t -4>xx- G (0) = F p (x, t) + F c (x) . (9) 

Equation (9) represents a system with explosive behavior as that shown in fig- 
ures 1,2 (a) and 3 (a), when the control perturbation F c (x) = 0. The problem 
is to find a controlling perturbation F c (x). Suppose F p (x, t) is a function that 
periodically generates the instability conditions discussed in the first part of 
the paper. An example can be the following F p (x,t) = a cos (cut) tanh (Bx). 
The strategy could be to find a perturbation F c (x) such that the superposition 
F p (x, t) + F c (x) does not satisfy the instability condition anymore for any t. 

It is remarkable that this can be achieved using a localized perturbation. For 
instance F p (x,t) = —0.385 cos (cut) tanh (Bx) is a turbulent-producing per- 
turbation, and F c (x) = 0.75 sinh (Bx) cosh -3 (Bx) can control this behavior. 
This can be seen in Fig. 3 (b). 

Similarly, the turbulence created by nonlinear damping can be controlled with 
a stabilizing perturbation: <p u — <Pxx — b(pt + o0 3 — G ((f)) — F c (x) . 

The explanation for these phenomena is based on our analytical results pre- 
sented above. The perturbation F p (x,t) = —0.385 cos (cut) tanh (Bx) is able 
to destroy the soliton and produce a highly nonstationary state because the 
perturbation F(x) = —0.385 tanh (Bx) leads to the instability of the soli- 
ton shape mode (according to our condition (5)). So F p (x,t) will produce 
this condition regularly. The soliton will be exposed to this instability again 
and again. The soliton destruction produces several new solitons and anti- 
solitons which also can be later destroyed because the perturbation makes 
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them unstable too. The control perturbation F c (x) is able to stabilize the 
soliton because the total perturbation F (x,t) = —0.385 cos (out) tanh (Bx) + 
0.75 sinh (fix) cosh~ 3 (fix) does not satisfy the shape mode instability condi- 
tion for any t. That is, when we solve the stability problem (3) for F(x) = 
— 0.385/xtanh (Bx) + 0.75 sinh (fix) cosh -3 (fix), the internal shape modes are 
always stable for — 1 < /j, < 1. 

The kink-solitons are examples of a very general phenomenon called topologi- 
cal defects. This set of phenomena includes: topological solitons, vortices and 
spirals [11]. Although these objects can possess different origin and nature in 
different physical systems, they all possess very similar dynamical properties 
[11]. 

The breakup of topological defects has been observed in experiments in many 
systems [12]. 

In different experiments, it has been observed that one topological defect 
can breakup into several topological defects. In particular, the "elementary" 
breakup that we have found, where one topological defect breaks up into three 
topological defects: one antidefect and two defects, has been observed in car- 
diac tissue [13]. 

All the situations discussed in the present paper that lead to very complex 
spatiotemporal behaviors start with soliton breakups (see figures 1, 2 (a), and 
3). 

At least, the following two different breakup scenarios are documented in 
experiments [14,15,16]. In one case, the breakup (leading to turbulence) occurs 
when a spatiotemporal external forcing is added to the system [14]. In a second 
case, the topological defects break after a Hopf bifurcation [16]. 

We believe that the results of the present paper show that very similar phenom- 
ena can occur with kink-solitons in Klein-Gordon systems. We have been able 
to produce defect-mediated turbulence using spatiotemporal external forcing, 
and after a Hopf bifurcation generated by nonlinear damping. 

A. Bellorin would like to thank OPSU-CNU for support. 
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